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MULTIMODAL FAR-FIELD ACOUSTIC RADIATION PATTERN - AN APPROXIMATE EQUATION 


Edward J. Rice 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 4i4135 


Abstract 


The far-field sound radiation theory for a 
circular duct was studied for both single mode and 
multimodal Inputs. Tlie investigation was Intended 
to develop a method to determine the acoustic power 
produced by turbofans as a function of mode cut-off 
ratio. This information is essential for the design 
of acoustic suppressors in engine ducts* With rea- 
sonable simplifying assumptions the single mode ra- 
diation pattern was shown to be reducible to a func- 
tion of mode cut-off ratio only (modal indices re- 
moved) . With modal cut-off ratio as the dominant 
variable, multimodal radiation patterns can be re- 
duced to a simple explicit expression. This approx- 
imate expression provides excellent agreement with 
an exact calculation of the sound radiation pattern 
using equal acoustic power per mode. Radiation pat- 
terns for cases other than equal modal power are 
presented using the approximate radiation equation. 
An approximate expression for the duct termination 
losses as a function of cut-off ratio is also in- 
cluded. 


Introduction 


The far field acoustic radiation pattern pro- 
duced by a single imde in a circular duct has been 
presented for flanged ductsl»2,3 as well as for un- 
flanged ducts. Multimodal radiation patterns 
representing the summation of the contributions of 
the many modes have been studied.^ Besides being 
very time consuming it was not obvious how cases 
other than equal acoustic power or amplitude per 
mode could be handled using exact calculations for 
the multimodal case. A method is presented here 
which makes this calculation possible by using a 
cut-off ratio biasing function. 


I This paper presents an analysis in which the 

single mode radiation equations can be summed ana- 
lytically into a multimodal radiation pattern. The 
single mode theory is first simplified into a func- 
tion only of mode cut-off ratio which means that 
^ the actual modal identification is suppressed, and 
i all modes with approximately the same cut-off ratio 

* will have approximately the same radiation pattern. 
A modal density function, which is also a function 

i only of mode cut-off ratio, is used to estimate the 
number of modes significantly contributing to the 
pressure at a far-field observation point. The 
single mode pressure amplitude and number of modes 
contributing are then combined to give a resultant 
! mean-square pressure at the far-field observation 
i point. This procedure yields a simple explicit ex- 
I pression which is then compared to exact calcula- 
; tlons using equal power per mode. The extension 

• to the more general case of unequal power per mode 
is also presented. An approximate expression for 
the duct termination loss as c* function of cut-off 
ratio is developed. Calculation.' using the multi- 
modal far-field radiation expressions are compared 
with experimental data. 

The concept that the far-field radiation pat- 
tern is a function of acoustic power versus mode 


cut-off ratio Is useful for acoustic liner design 
purposes. The propagation of sound through lined 
ducts has previously been shown to be a function of 
cut-off ratio. 7 The work presented here repre- 
sents a part of a simplified acoustic liner design 
program under development. The noise source infor- 
mation needed for the liner design can be obtained 
directly from the far field radiation pattern. 
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mic 


speed of sound, m/sec 
duct diameter, m 

modal density function, see eq, (15) 
frequency , Uz 

acoustic intensity tinnsmitted out from 
termination, N/msec 

acoustic intensity incident upon duct ter- 
mination, N/msec 

derivative of the Bessel function of the 
first kind of order m with respect to 
its argument 

side lobe index (see eq. (24)) 

spinning mode lobe number (also order of 
Bessel function) 

fraction of the total number of modes whose 
principal lobes contribute at any far 
field position, see eq. (16), 

exponent on cut-off ratio biasing funOftion, 
see eq. (18) 

far field acoustic pressure amplitude, see 
eq. (1), N/m^ 

acoustic pressure in the duct, N/m^ 

P due to multimodal summation of principal 
lobes, N/m^ 

r at the peak of the principal lobe, N/m^ 

2 

P due to side lobes, N/m 
reflection coefficient at duct termination 
radial coordinate in cylindrical system, m 
radius of cylindrical duct, m 
time, sec 

hardwall duct mode eigenvalue, see eq. (3) 
Cut-off ratio biasing function, see eq. (18) 
frequency parameter, fD/c 

duct termination dimensionless acoustic 
radiation resistance 

radial mode number 

mode cut-off ratio, see eq. (4) 

cut-off ratio for mode whose principal lobe 
peaks at a given far field microphore 
position 
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T 
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cut-off ratio for a particular sidolobe 
peaking at a given angle, see eqs. (26) 
and (27) 

propagation coefficient, see eq^ (21) 

angular coordinate, radians 

duct termination dimensionless acoustic 
radiation reactance 

far field angle measured from inlet duct 
centerline, degrees 

^ for the peak of the principal lobe, see 
eq. (9) 

for the peak of the side lobes, see 
eq. (24) 

circular frequency, radians /sec 
Analysis Technique and Results 


The single mode far- field acoustic radiation 
expression will first be developed as a function of 
cut-off ratio. The summation required for the 
multimode case will then be derived. The equations 
without steady flow will be used here as an illus- 
tration. 


Single Mode Equation 


A far-field mean-squared pressure radiation ex- 
pression for a flanged circular duct for equal 
acoustic power per raodeS can be expressed as, 


(rn)^ |/(iTn) ^ -~ct^ sin^4< [J^ (irn sin 




( 1 ) 


[a^ - (ttd sin 


where the frequency parameter q Is, 

fD 


( 2 ) 


and the mode eigenvalue ot is determined from 


j;(a) ^ 0 


(3) 


Equation 1 was derived for zero tlach number, negli- 
gible end reflections, and a constant has been 
omitted. Note that subscribts m,p should occur 
on and a to denote the lobe number m and 

the radial mode index p, but they are omitted here 
for brevity. 

Let the mode cut-off ratio be defined as, 
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(4) 


which is consistent with previous results. 2,8 Using 
eq. (4) in eq, (1) yields, 
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m 

[l/s^ - sin^4»l^ 


(5) 


Note that in eq, (5) the modal identity remains in 
the m/a term and in the Bessel function order m. 
Two assumptions are made. First, 



( 6 ) 


and second. 


J» (irn fiin 4<) ^ ^ Cirn sin tfi - a) (7) 

® Y TT^n sin ^ 


which can be obtained from published expressions.® 

Equations (6) and (7) are valid except for the 
first few radial orders of the higher lobe number 
modes- Using eqs. (4), (6), and (7) in eq. (5) 
yields , 

2 ^ 2 sin - l/g^{sln [Trn(sin if) - l/g)]}^ 

TT^n [1/C^ - 3±n\]^ 

( 8 ) 


Equation (8) shows that the radiation pattern is a 
function of cut-off ratio (?) only, and that modal 
identity has been suppressed being now implicitly 
involved in The frequency parameter (q) of 

course remains but this Involves only the usual in- 
puts of frequency and duct diameter (D) . Provided 
that conditions given by eqs. (6) and (7) are not 
radically violated, all modes with the same cut-off 
ratio have identical radiation patterns. One dif- 
ference between the approximate eq. (8) and the more 
exact eq. (1) occurs in the number of side lobes 
occurring forward (^ < ) of the principal lobe. 

Equation (8) shows that ror a given cut-off ratio 
(^) and frequency parameter (q) there are a constant 
number of these side lobes. Howevi^r, for eq. (1) 
the number of these side lobes is dependent upon 
the radial mode number (y) . This error occurs due 
to the difference between the small argument be- 
havior of the sine function and the Bessel function. 
Since the side lobes are generally ten to fifteen 
decibels below the principal lobe this error ia the 
forward side lobes is of little conoecaence in a 
multimodal radiation pattern. 

From eq. (8) it can be seen that both the nu- 
merator and the denominator go to zero when 


sin 


1 
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(9) 


which approximately defines the location of the 
principal lobe peak pressure. The subscript (p) has 
been added to to signify the angle for the peak 
pressure. A limiting procedure on eq. (8) yields 
for the magnitude of the principle lobe, 

p2 ^ ne V . ^ u sl = (10) 

P 2 2 sin .|.p 

where eq, (9) was also used. The same result occurs 
using the exact eq. (1) if the approximation of 
eq. (6) is used. 

Equation (8) also shows that the principal lobe 
extends between angles 
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sin 1(1 it I - i (11) 

and 

aln i|j % ■- + - (12) 

t n 

which represent the first zeros of 
sin [rnCsin ^ - l/5)l 

about the peak. 

Multimodal Radiation Pattern 
Equal Acoustie Power Per Mode 


which represents the fraction of the total number 
of modes whose principal lobes contribute to the 
pressure at angle The simplest possible way of 
integrating the amplitude and mode number effects 
is to multiply eq, (10) by eq. (16) to yield 

% 2 COB V- (17) 

Equation (17) is scon to be a simple exprcaalon 
for the far-fiold pressure at any angle due to 
many modes with equal acoustic power per mode. In 
figure 2, cq. (17) is compared to an exact calcula- 
tion using an exact auraniatlon of eq. (1) for over 
one thousand modes (A. V. Saule» unpublished), the 
agreement between the exact and approximate calcul- 
ations is seen to be excellent. Tills justifies the 
approximations made earlier. 

Unequal Acoustic Power Per Mode 


Consider a far-ficld location, say a microphone 
location, described by the angle With a multi- 

tude of modes present there is some mode with o cut- 
off ratio ^ whose principal lobe peaks at this 
location* The principal lobe of this mode is shown 
as the solid line with cut-off ratio shown us Kmlz 
in figure 1. Also there are many other modes whose 
peaks are not at ^[^ but whose principal lobe still 
contributes to the pressure at i|/* Samples of such 
modes arc shown as the dashed lines in figure 1. 

Tlie dash-dot curves show the extreme values of cut- 
off ratio beyond which no contribution to the micro- 
phone pressure la made by the principal lobes. The 
range of cut-off ratios which have some principal 
lobe contribution at i|j can be found from eqs. (11) 
and (12) to be 

i sin - - (13) 

to 

% siu ||« + - (1^) 

h 


It is desirable to have a simple oxpresaion 
comparable to equation (17) for other than equal 
acoustic power per mode* It is also of interest to 
alter the expression so as to keep it a function of 
cut-off ratio since propagation In acoustic liners 
has been correlated with cut-off rQtlo.^»7 A cut- 
off rotio biasing function can be defined as 


l/c” 


(IB) 


where the exponent (n) sexrves to alter the acoustic 
power distribution as a function of cut-off ratio 
which represents a modal power shift. Equation (1) 
can be multiplied by the cut-off ratio biasing func- 
tion l/5®^ to produce a shift In acoustic power 
away from equal power pev mode (n » 0). This is not 
the only biasing function which could be used, but 
it is both simple and convenient in form niid will 
be used in this initial effort. Exact colculatlons 
of tlic sum of the me&n-square pressure at any micro- 
phone position (fixed Vi) could be made by summing 
over all the propagating modes (varying eigenvalue 
a) with the cut-off ratio used for each mode to cal- 
culate the power biasing function. 


The number of modes in this cut-off range can be 
estimated using the modal density function, 

(15) 

r 


As in the case of equal acoustic power per 
mode, instead of performing the exact calculation 
Just outlined, a simpler procedure will bo used to 
approximate the summation over modes or Integration 
over cut-off ratio. First eq. (10) fqt the princi- 
pal mode peak value is multiplied by the cut-off 
ratio biasing function to yield. 


Tills is a density function in the usual sense in 
that an integration must be performed to obtain the 
number of modes. The most desirable way of obtain- 
ing the microphone pressure signal Is to multiply 
eq. (8) by the modal density function and integrate 
Over cut-off ratio between zero and infinity. Since 
this cannot be done in closed form a simpler proce- 
dure will be used. Integration of eq. (15) between 
the limits given by eqs, (13) and (14) yields 


^ - 1/5^ 09) 

P 2 27 


Equation (19) is then multiplied by the number of 
modes (eq. (16)) whono principal lobes can contri- 
bute to the mean-square pressure at the far-fleld 
angle iju This results in, 



4 ^ 4 sin lit 


(16) 


P** 2 cos iji sin 

m 


X20) 


Figure 3 shows a plot of eq. (20) for fouv values 
of the exponent n. As n is increased the acous- 
tic powor distribution is shifted coward cut-off 


(5 “ 1) and the radiation pattern peaks toward 
higher angles with a noticeable reduction in the 
power near the inlet axis. The peak value of 
eq. (20) is located at tan 

The biasing function (l/s”) used here is de- 
sirable for n ^ 0 but is not completely satlsfac- 
tory for n < 0* In the latter case no absolute 
maximum occurs and the mean-square pressure near the 
Inlet axis continually increases as n is reduced 
(mode larger negative) . As mentioned previously 
this first attempt at the biasing function was made 
for convenience and further experimentation must be 
performed when making quantltiitive fits to real 
data. Qualitative fits to several s^ts of data will 
be shown after reflection from the duct termination 
effects and the contribution of side lobes have been 
Included in the radiation pattern. 

Duct Termination Loss as a Function 
of Cut-Off Ratio 


The acoustic power loss (due to reflection) at 
the duct termination will be determined for each 
individual mode. Thus modal interactions as con- 
sidered by Zorumski^^ will be neglected and the re- 
sults are only approximate. The power loss will be 
shown to be a function only of cut-off ratio. 


The pressure in the circular duct has the form, 
-iuiTx +imTx1 


Pd - e 


imt-imi^) 


+ Re 


( 21 ) 


where 


ations among modes in the oscillatory parts of the 
curves beyond the cut-off point. This implies that 
the radiation impedance is approximately dependent 
only on cut-off ratio. Also note in eq. (23) that 
the remaining term (t) is only a function of cut- 
off ratio (see eq. (22)), Thus the duct termination 
loss is approximately only a function of cut-off 
ratio. 

Calculations using eq. (23) and the radiation 
impedance components of figure 4 are shown in fig- 
ure 5. The slight difference in radiation impedance 
components causes a small difference in termination 
loss, A further approximation will be made com- 
pletely ignoring the small differences in radiation 
impedance among modes (when considered as a function 
of cut-off ratio;. For relatively large cut-off 
ratios 1 and Xr 0* Using these values in 

eq. (23) yields. 


m) 

INC (1 + t) 


Equation (24) is purely a function of cut-off ratio 
and is plotted versus (^ - 1) as the solid line in 
figure 5. This approximate relationship Is seen to 
adequately describe the termination loss calcula- 
tions shown in figure 5 and due to the previously 
described results of reference 12 is expected to 
provide an adequate result for all of the modes. It 
can be seen in figure 5 that the duct termination 
loss is large only near mode cut-off. For cut-off 
ratios greater than 1.2 the loss is less than one 
decibel. Equation (24) has been used for the ter- 
mination loss for the results which are presented 
in later figures. 




given in terms of mode eigenvalue and cut-off ratio. 
Note that hard walls and zero axial Mach number are 
used here for simplicity. The reflection coeffi- 
cient R in eq. (21) accounts for the reflection 
back from the duct termination which represents 
acoustic power lost to the far-field. The boundary 
condition nt the duct end, in terms of radiation 
resistance Op and reactance Xr> used to eval- 
uate R, The ratio of transmitted to Incident 
acoustic power Intensity can then be evaluated as 


llNC (1 + + (tx^)^ 


The radiation impedance components (0^ and Xr) 
the first two axisymmetric (m *» 0) radial modes are 
plotti ^ figure 4, These calculations were made 
usj^W^ method of reference 11- Note that when 
calculated impedance components are plotted 
Versus mode cut— of f ratio the results for the two 
Uiodes are aXcCSt coincident. In addition, radiation 
Impedance calculations for spinning modes were pre- 
sented in reference 12 plotted as a function of 
nondiiticu^ianal frequency minus cut-off frequency* 
Thes0 calculations shnVJed the same trends as shown 
in tiRUlC Ax n^^tdance components rapidly In- 

xit%hm (C Vith only minor vari- 


The termination loss ^hus joins the list of 
duct lining phenomena which appear to be approxi- 
mately related only to mode cut-off ratio. These 
Include acoustic liner modal optimum impedance and 
maximum possible attenuation^*' and modal far-field 
radiation pattern (eq. (8)). 

Effect of Sidelobes on Multimodal 
Radiation Pattern 

As previously discussed, the sidelobes In the 
radiation pattern have considerably lower amplitude 
than the principal lobe and were ignored in the 
approximate calculations presented in figures 1 and 
2. They were, of course, automatically contained in 
the exact calculations of figure 1 and are seen to 
have very little effect. However, when the termina- 
tion loss is included which mainly reduces the 
principal lobes at higher angles (near cut-off), the 
radiation pattern is dominated by the sidelobes from 
modes whose principal lobes fall nearer to the inlet 
axis (well propagating modes with larger cut-off 
ratios). Also when duct liners are included, the 
principal lobes may be drastically reduced over a 
far-field angle range, and a noise floor results 
from the sidelobes of the less attenuated modes. 

The sidelobes are thus included in the results shown 
later in this paper. 

Equation (8) shows that the k^^ sldclobe 
(counting from the principal lobe) peaks at the 
angle (li^pg) given by, 
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I 


sin ** ^*^*^* * • 

wherg thci plus and minus signs take earn of side- 
lobes at angles greater than and leas than that of 
the principal lobe. Equation (?.5) is the parallel 
of eq. (9) for the principal lobe. Using eq. (25) 
and using a procedure similar to that used for the 
principal lobe the number of modes vhoac side- 
lobe contributes to the mean-square pressure at the 
far field angle ^ is determined as^ 


Figure 7 shows some broadband datal^ for the 
NASA Fan B on o static outdoor test stand. This 
data has been carefully analysed to remove the 
multiple pure tones from the broadband contribu- 
tion.^^ Also an attempt has been made to remove the 
radiated noise contributed from the fan aft duct or 
jeti A comparison of the data with the equal acous- 
tic power per mode calculation shows the data to a 
agree very well with this curve. Thus for the two 
different fans tested in two different facilities, 
as represented by figures 6 and 7, the broadband 
noise radiation patterns appear to be very close to 
that expected for equal acoustic power per mode. 


N 


S 



(26) 


where 5ps is the cut-off ratio of the mode whose 
k^h fiidelobe peaks at the angle and is given by 


= sin k ^ n sin - 1/2 (27) 


for sidelobes above the principal lobe and 


P® (28) 

for sidelobes below the principal lohe. 

Finally the multimodal contribution of the k*^^ side- 
lobe is 


A sin » 

” 2 2,n+l r, ,,2 . 2 .12 


(29) 


A summation over the sidelobes is required using 
eqs. (27) and (28). 


Blade passage frequency data from two fans 
statically flt two speeds ara shown in fig- 

ure 8, Although there are some scatter in the data, 
the trends fit relatively well with the theoretical 
curves shown with n “ 1/2 or 1. These data are 
shifted away from equal power (n ® 0) towards cut- 
off. This is evident from the higher radiation at 
larger angleB and the reduced radiation toward the 
inlet axis. 

A i) at pie of the blade passage frequency di- 
rectivity data "" r the Lycoming YF-102 engine tested 
at NASA-Lewis is shewn in figure 9. The solid curve 
is the multimodal theory for equal acoustic power 
per mode. Except for small aberrations near the 
axis and centered at 65°, the curve fits the data 
very well. The devlatitin at 65^ may represent the 
contribution of an additional single mode, the rotor 
locked mode wlt^: forty lobes (m = 40) . Only the 
first radial order of this mode can propagate. The 
dashed linjC represents the adding of this mode to 
the theoretical curve. This helps to explain the 
deviation the higher angles. This illustrates 
a reasonable procedure which can probably be used 
in future data analysis. After deviations in the 
radiation pattern have been determined to be truly 
valid, a limited number of physically sensible modes 
may be added to the underlying smooth multimodal 
radiation pattern. No attempt has been made to ex- 
plain the deviations of the data near the axis. 


Comparison of Calculated Radiation 
Patterns with Experimental Data 

Several sets of experimental radiation direc- 
tivity patterns will be compared to the calculated 
patterns derived in the previous sections. Although 
the ultimate intent of this work is to infer the 
acoustic power distribution as a function of cut-off 
ratio, this has not been done rigorously here* The 
cut-off ratio biasing function is considered pre- 
liminary and for Illustrative purposes only at this 
time. Several interesting qualitative features of 
the data-theory comparisons, however, will be 
pointed out. 

In figure 6 four sets of data considered to be 
broadband nolsei^ are plotted. Only inlet noise is 
present since aft noise is ducted away. The solid 
curve represents the multimodal equal power per 
mode theory with the biasing exponent equal to zero 
(n =* 0) « The comparison of the data with the theory 
shows this broadband noise to be very close to equal 
acoustic power per mode. The data arc seen to be 
almost unchanged when the tunnel flow Is added. 

Equal power per mode is a very appealing result for 
broadband noise since this is usually associated 
with random processes. 


This very limited comparison of data with 
theory has pointed out that a variety of acoustic 
power versus cut-off ratio (or modal) distributions 
are to be expected from experimental data. Some of 
the dlffereneces might be due to different fan de- 
signs or to different experimental test installa- 
tions . 

ConcludlnR Remarks 

A simple explicit equation has been derived for 
multimodal far field acoustic radiation from a 
flanged circular duct for arbitrary acoustic power- 
cutoff ratio distribution. The approximate expres- 
sion agrees very well with exact calculations when 
equal acoustic power per mode is specified. The 
approximate expression was based upon only the 
principal lobes of the modal radiation pattern, but 
additional expressions arc provided to account for 
the sidelobes in cases where they must be included. 
The derivation used an approximate integration of a 
continuously distributed modal distribution which 
will certainly be of questionable validity if only 
a few modes are propagating. 

Approximate expressions for the single mode 
far-field radiation pattern and duct termination 
loss were presented. Exact modal identification 


5 


hos been suppressed in these two expressions and 
they have been shown to be functions only of the 
icodal cut-'off ratio. This is convenient since pre- 
viously it has been shown that the sound propagation 
in the noise suppressor itself is predominantly only 
a function of cut-off ratio. 

Several sets of experimental data for inlet 
radiated noise were compared to theoretical curves* 
This comparison^ although qualitative at this time, 
has shown that a variety of acoustic power-cutoff 
ratio distributions might be expected. These may 
depend upon the spectral nature (tone or broadband) 
of the noise, the fan characteristics, or the test 
Installation. Some distributions are fairly smooth 
versus angle while some may require the addition of 
some physically sensible modes. 
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Figure 1. - Schematic of modal principle lobes of radiation 
contributing to a far*field microphone signal. 
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EXACT CALCULATIONS USING 
BESSEL FUNCTIONS 

APPROXIMATE CALCULATIONS 



Figure 2 . - Comparison of exact and approximate multi- 
modal far-field directivity patterns for equal acoustic 
power per mode. 



o □ 


Z UJ 
LU ^ 
Q- 3 
O ti- 


□ 


o □ 


. oo 

- 'O 


^ J_i_J I L 


.-4 OO O KT 


ao ^ ‘33NV13V3a 

ao 33Nvisisia Noiiviava NoiivNiwa3i i3na 



sp 'i3/\ri 3anss3ad ONnos 3Aiitn3a 


ANGlf FROM INLET AXIS, deg CUT-OFF RATIO. E 

Figure 3. - Far field multimodal directivity variation with Figure 4. - Variation of duct termination radiation imped- 

cut-off ratio bias function exponent. ance components with cut-off ratio. 
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Figure 3. - Modal transmission loss at duct termination. 
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Figure 6. - Comparison of multimodal equal power per 
mode directivity to experimental data obtained in a 
wind ti'nnel (ref. 13). 
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Figure 7. - Compctrison of equal power per mode multi- 
modal directivity to experimental data from NASA 
Fan B (ref. 14). 
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Figure 8. - Comparison of fan blade passage frequency 
data with calculated far-field multimodal directivity 
patterns. 
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Figure 9. - Comparison of experimental data with cal- 
culated multimodal directivity pattern. 
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